Universality of Fluctuation-Dissipation Ratios: The Ferromagnetic Model 
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We calculate analytically the fluctuation-dissipation ratio (FDR) for Ising ferromagnets quenched 
to criticality, both for the long-range model and its short-range analogue in the limit of large 
dimension. Our exact solution shows that, for both models, X°° = 1/2 if the system is unmagnetized 
while X°° = 4/5 if the initial magnetization is non-zero. This indicates that two different classes of 
critical coarsening dynamics need to be distinguished depending on the initial conditions, each with 
its own nontrivial FDR. We also analyze the dependence of the FDR on whether local and global 
observables are used. These results clarify how a proper local FDR (and the corresponding effective 
temperature) should be defined in long-range models in order to avoid spurious inconsistencies and 
maintain the expected correspondence between local and global results; global observables turn out 
to be far more robust tools for detecting non-equilibrium FDRs. 



I. INTRODUCTION 

One of the main goals of modern statistical mechanics 
is to find a general theory of non-equilibrium processes. 
Although significant advances have been made in the 
past [1, 2], a complete theory of non-equilibrium systems 
analogous to thermodynamics does not yet exist. A com- 
mon approach has been to extend well-known equilibrium 
concepts to the non-equilibrium regime. One of the most 
important among these is the temperature, and many re- 
searchers have therefore tried to extend it to systems out 
of equilibrium by introducing a so-called effective temper- 
ature [3] . The question of whether and how such a quan- 
tity can be defined properly is a central issue in the con- 
struction of a general theory of non-equilibrium systems. 
Different definitions have been employed in granular sys- 
tems [4] , driven systems [5] and glassy systems [6] among 
others. Studies of mean-field spin-glasses have shown 
that an effective temperature can be defined by measur- 
ing the violation of the fluctuation-dissipation theorem 
(FDT) in terms of a so-called fluctuation-dissipation ra- 
tio (FDR) [7] 



X{t,t^) 



TR{t,t^) 

dC(t,t„) 



where T is the temperature of the heat bath, C{t,t^) 
is the autocorrelation function of a given observable and 
R(t,t^) the conjugate response function; the latter en- 
codes the change of the value of the observable at time t 
to a small perturbation at an earlier time t„. In equilib- 
rium FDT is verified and X{t, iw) = 1- Whether the FDR 
is useful more generally, and in particular beyond mean- 
field models, has been the subject of debate in recent 
years [8, 9]. R has been shown that the effective temper- 
ature defined through the FDR, Tcff (t, iw) = T/X{t, t^), 
has good thermodynamic properties for some mean-field 
models [6]; the zeroth law of thermodynamics can also 
be extended to the non-equilibrium regime [10]. Never- 



theless, there are still many open questions regarding the 
physical meaning and universality of the FDR [11, 12]. 

An important aspect in the study of FDT violation 
in glassy systems is its application to ferromagnetic sys- 
tems which are quenched from high temperature to the 
critical temperature (see e.g. Refs. [13-15] and the re- 
cent review [16]) or below. The ensuing non-equilibrium 
evolution, where the system coarsens - by the growth of 
domains with the equilibrium magnetization, for T < Tc 
- is of course different from that of glasses in many re- 
spects; for example, thermal activation effects are irrele- 
vant for the long-time dynamics. However, there are also 
appealing similarities. In particular, equilibrium is never 
reached in an infinite system and as a consequence the 
system exhibits aging, i.e. a dependence of the relaxation 
properties on the time elapsed since the quench. 

The simplest ferromagnetic model that can be studied 
at criticality is the Ising chain first solved by Glauber [17]. 
At the critical point, T = 0, the magnetization jumps dis- 
continuously from (T > 0) to 1. The relaxation dynam- 
ics at zero temperature after a quench from T = cxd, i.e. 
a random initial configuration, has been studied in e.g. 
Refs. [18-20]. Recently, the FDR has been calculated an- 
alytically for a randomly staggered perturbation and the 
corresponding spin autocorrelation function [21, 22]. For 
long times t and tw this gives X{t,t^) (1 + iw/^)/2. 
In the limit t ^ oo the FDR then approaches X — 1/2, 
which coincides with the value obtained in models char- 
acterized by diffusive dynamics (such as the random walk 
or the Gaussian model [23]). These results have lead to 
the suggestion [13, 24] that for systems at criticality the 
limiting value of the FDR 



= lim lim X{t, t^ 

— ^OO t — 'OC 



is a universal quantity. Consistent with this, the exact 
solution of the ferromagnetic spherical model in d di- 
mensions at criticality also gives X°" — 1/2 ford > 4, 
i.e. above the upper critical dimension [13, 24]. 
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In considering the universality of the hmiting FDR, one 
issue is whether and how the limiting FDR depends on 
the observable whose correlation and response are mea- 
sured. An obvious alternative to the local spin autocor- 
relation is its long-wavelength analogue, i.e. the corre- 
lation function of the fluctuating magnetization. Exact 
calculations for the Ising chain [14, 25] and the spheri- 
cal model [26] as well as numerical simulations [14, 27] 
for the Ising model in = 2 show that the resulting 
global X°° is always identical to the local version. This 
local-global correspondence, which can also be obtained 
by field-theoretic arguments [15, 16], is rather reassuring; 
physically, it arises because the long wavelength Fourier 
components of the spins are slowest to relax and dom- 
inate the long-time behaviour of both local and global 
quantities. For a numerical determination of the limiting 
FDR the global quantities are often more suitable [14, 27] 
because plots of susceptibility (integrated response) ver- 
sus correlation are close to straight lines with slope X°°. 
We do not discuss in this paper the FDR for other ob- 
servables that are nonlinear in the spins, e.g. the energy; 
this question is studied in [15, 16, 26]. 

A further key question is to what extent the limit- 
ing FDR X°° depends on the initial conditions for the 
coarsening dynamics. The results quoted above all ap- 
ply to initial high-temperature equilibrium, i.e. an un- 
magnetized system with no or only short-range spatial 
correlations. Some initial progress in considering more 
general initial states has alrc^ady hnm made. For the 
Ising chain [28] a non-zero initial magnetization does not 
change the value of X°°] this is unlikely to be a general 
result, however, because the Ising chain at its T = crit- 
ical point has the peculiarity that the magnetization re- 
mains constant instead of decaying to zero. Other values 
of the limiting FDR are nevertheless possible even in the 
Ising chain; e.g. when the magnetization is zero initially 
but correlations between spins are so strong that only a 
finite number of domain walls exist in the system, one 
finds [29] X°° = 0. A more general analysis for unmag- 
netized but spatially correlated initial conditions in the 
spherical model gives qualitatively similar results [28]: if 
initial correlations are strong (i.e. decay with distance as 
slow power laws), the limiting FDR is =0; otherwise 
it is the same as for an uncorrelated initial state. 

In this paper wc investigate in full detail the depen- 
dence of the FDR on the initial condition for a ferromag- 
net with long-range interactions as well as short-range 
interactions in the limit of large dimension d. Our main 
result will be that the non-equilibrium dynamics starting 
from unmagnetized and magnetized initial conditions are 
in different universality classes that are distinguished by 
different and nontrivial, i.e. non-zero, values oi X°° . We 
also analyze the correspondence between local and global 
FDRs. In a naive analysis this appears broken in the 
long-range case, but we show that it is recovered when 
finite-size corrections are included. 

The paper is organized as follows. Section II describes 
in outline the calculation of the basic evolution equations 



for correlation and response in the long-range ferromag- 
net. In section III wc compute from these the global 
and local FDRs. Section IV contains the correspond- 
ing analysis for the short-range ferromagnet, where the 
lengthscale-dependence of the FDR can be made explicit 
and the local-global correspondence holds as expected. 
In section V this correspondence is shown to hold also 
for the long-range model once finite-size corrections are 
accounted for. Section VI, finally, summarizes our key 
results and conclusions. Technical details are relegated 
to two appendices. 



II. LONG-RANGE FERROMAGNET 

In this section we study the ferromagnet with long- 
range interactions and in particular the non-equilibrium 
dynamics of the relaxation functions (correlations and 
responses). We leave technical details to Appendix A and 
focus here on the conceptual aspects of the calculation 
and the results. 

The model is defined by the Hamiltonian 

where the ai = ±1, i — \ . . .N , denote Ising spin vari- 
ables and hf^^ is a position-dependent external field. The 
strength of the coupling between spins, Jo, is normalized 
by a factor 1/{N — 1) to ensure an extensive energy; we 
will take Jq = 1 without loss of generality. The dynamics 
we consider is of Glauber type: each spin Ui flips inde- 
pendently with rate [1 — crj tanh(/3/ij)]/2 where (3 = 1/T 
is the inverse temperature and hi the local field acting 
on spin i, 

^^=''^' + N^ll2''^- (2) 

To keep the notation compact we define the function 
th(z) = ta.nh{(3z) and abbreviate ti = tanh(/3ft,j) = 
th(/ii). Multiplying the flip rate by the change —2ai when 
(Ti flips, it follows that 

d 

Q^{(^i) = {ti - (Ti) (3) 

where the brackets denote an average over the thermal 
history of the system and over the initial conditions. For 
spin products {i ^ j) one finds similarly 

d 

Q^icTiaj) = {{ti - (Ji)(jj) + {ai{tj - Gj)). (4) 

Throughout this paper, omitted time arguments indicate 
dynamical averages evaluated at time t. To obtain the 

two-time correlation functions we can take advantage of 
the fact that (3) does not depend on the initial condition 
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and so is equally valid for correlations with some quantity 
at an earlier time < t. This gives 



d_ 

dt 



{(Ti{t)aj{t^)) = {[ti{t) - ai{t)]aj{t^)). (5) 



Using (3) we can also express the time evolution of the 
magnetization, m = AT"^ ^ 



dm 1 \ ^ , , 



(6) 



The relaxation dynamics of physical systems is routinely 

characterized by correlation and response functions. We 
will consider here the properties of both global and lo- 
cal relaxation functions. Our main goal is to clarify the 
relationship between them and the implications for the 
corresponding FDRs. Note that while the fluctuations of 
individual spins, as encoded in the local (auto-) correla- 
tion function, are 0(1), those of the fluctuating magne- 
tization A'^"^ 'J2i Ci are 0{N~^^^). Nevertheless we will 
see that the same physics can be extracted from both 
quantities. The global correlation function is scaled by 
a factor of A'' below to give an 0(1) quantity as in the 
local case. 

Previous studies [30, 31] have shown that the relevant 
correlation functions are the connected ones, defined by 



First we analyze the equal-time correlations. For i ^ j it 
follows from (4) that 



d_ 

at 



Cij{t,t) = {{t 



ai)aj) + {(Ti{tj - aj)) 

- {U-ai){aj) - {a^){tj -a^) 

2C,j(t, t) + (At.AcTj) + {Aa.Atj) (7) 



where we use the notation Aijj = — for the deviation 
of any quantity from its average, so that e.g. Ati = ti — 
{ti) and similarly for AcTi. For the two-time correlations 
one gets similarly from (5): 



d_ 
dt 



Cij{t,t^) = {[ti{t) - ai{t)]aj{U)) 

= -Cij{t,t^) + {Ati{t)Aaj{t^)). (8) 



Equations (7) and (8) are general and also valid for short- 
range systems provided that the appropriate local field 
replaces the long-range expression (2). 

To make progress, we exploit the fact that in the long- 
range model, for large iV, the correlations between differ- 
ent spins are of 0{N~^). So the fluctuations Ahi of hi 
aroimd its mean hf^^+m are small, of 0{N~^/'^). In equi- 
librium away from criticality one can indeed show that 
{{Ahif) = 0{N-'^) while {{Ahi)'^Aaj) and {{Ahif) are 
0{N~'^). We will assume that the correlations out of 
equilibrium are of the same order. This is reasonable if 



we start from an initial state with weak correlations and, 
for quenches to criticality, also restrict ourselves to the 
interesting non-equilibrium regime where all times are 
short compared to the equilibration time. 

Setting the external field to zero, we can Taylor expand 
the nonlinear terms in our equations of motion in powers 
of Ahi-. 

U = th(m) -I- Ahith'im) + ^{Ahifth"{m) + ... (9) 

Since we are only interested in the leading terms, we trun- 
cate this expansion after the linear term in Ahi; sublead- 
ing corrections are discussed in Sec. V. Since {Ahi) = 0, 
the leading order term in the equation of motion for the 
magnetization (6) takes the expected mean-field form 



dm , , . 

-— - = -m -I- th(mj. 
at 



(10) 



Due to spatial translation invariance (or more precisely 
permutation invariance) between spins in the long-range 
ferromagnet, there are only two different correlation func- 
tions, one local and one non-local. We write these as 



Cii = Cioc + 0{N-'') 



C\j=Cni/N + OiN-^) 



(11) 



(12) 



using the fact that the non-local correlations are only 
0(1/ N) to leading order. 

In terms of the local (11) and non-local (12) correla- 
tions, the global correlation is defined by 



= Cloc + (^nl 



(13) 



and gives the leading contribution of the correlator of the 
magnetization; see eq. (Al) of Appendix A. In order to 
compute the dynamical equation for the global correla- 
tion function (13) we need the equations for the local and 
non-local correlation functions, which can be expressed 
as 

-C„i(t,f) = -2aC„i(t,t) +2th'(m)(l -m^) (14) 



^Cni(i, iw) = -Cniit, iw) + th'(m)Cg(t, iw) (15) 



Cioc{t,t) = l-m'^{t) 



d_ 
dt 



C'loc(i)^w) — ~C'loc(i)^w) 



(16) 



(17) 



as shown in Appendix A (see eqs. (A2), (A5), (A6) 
and (A7)). The quantity a appearing in (14) is defined 
as 

a = 1 - th'(m) = 1 - /3[1 - ta.nh^{Pm)]. (18) 
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We will not normally write its time-dependence explicitly. 

In order to complete the analysis of the dynamics of our 
model, we need to find the linear response to applied ex- 
ternal fields. This is characterized by the response func- 
tions 



As for the correlation functions, we have to leading order 
in N 



Rii = Rioc + 0{N-^) 
Rij = Rni/N + 0{N-^) 



(19) 
(20) 



while the leading term in the global response, of the mag- 
netization to a uniform field, is 



-Rg — i?loc + Rnl- 



(21) 



The evolution equations for these response functions can 

be expressed as 



d_ 
dt 



Rnl{t, iw) = th'(m)[i?loc(t, tw) + Rnl{t, ^w)] " Rnl^t, 

(22) 

(23) 



d_ 

di 



Rlocii^ ^w) — Rloci,^: ^w) 



as derived in Appendix A (see eqs. (AlO) and (All)). 
These equations can be integrated forward in time start- 
ing from the values of the equal-time response func- 
tions (A13): 

Riocit, t) = th'(m), i?ni(i, t) = 0. (24) 

The evolution equations (14-17) for the correlations and 
(22-24) for the responses contain all the relevant informa- 
tion about the dynamical properties of the ferromagnetic 
systems. 



III. FLUCTUATION-DISSIPATION RATIOS 

Using the results obtained in the previous section, we 
can study in full detail the fluctuation-dissipation ratios 
(FDR) for global and local relaxation functions. Cru- 
cially, we will be able to investigate how the value of the 
asymptotic global FDR depends on the initial condition. 
Our calculation in the infinite system size limit will pro- 
duce different values for global and local FDRs; this ap- 
parent breaking of the expected local-global correspon- 
dence will be solved and carefully explained in Section 
V. 



A. Global FDR 

To compute the global FDR, we first need the equal- 
time global correlation Cg{t,t) = C\oc{t,t) + Cni{t,t) = 



1 — m^{t) -f Cni(t, t). Differentiating this expression with 
respect to time and using (10) and (14) we get 



d_ 

dt 



Cg{t,t) = -2m[th(m) - m] - 2aCni{t,t) 

+ 2(l-a)(l-m^) 
= -2aCg{t,t) + b (25) 



with & = 2[1 — mth(TO)] and a defined in (18). Equa- 
tion (25) can be integrated explicitly to get 

CS,t) = r'it)C,{0,0) + J'dt'^bit') (26) 



where we have defined the quantity 



r{t) = exp 



dt' a{t') 



(27) 



For the global two-time correlation (13), the sum of (15) 
and (17) gives {d/dt)Cg{t,ty,) = —aCg{t,tv,) and after 
integration 



Cg{t, ttfi) 



r{t) 
r(tw) 



(28) 



To compute the corresponding global response (21), we 
add equations (22) and (23) to obtain 

d 

— Rg{t,t^) = —aRg{t,tvf). 
The solution for the global response is then given by 

Rsit, *w) = ^r^-^tl - tanh'(/3m(tw))] (29) 

where we have used Rg{t,t) = Rioc{t,t) and (24). 

By combining the last two results we get an exact an- 
alytical expression for the global FDR 



Xg{t,ty,) 



T[l - a{t^ 



b{t^) - a(fw)Cg(iw,iw) 



(30) 



where we have also used the fact that dr{tyj)/dtyj — 
—a{tyf)r{t^). Equation (30) shows explicitly that the 
global FDR depends only on the earlier time t^, which is 
a feature often seen in simple mean-field models. From 
the general expression (30) we can now analyze the 
asymptotic FDR for different initial conditions. 



1. Zero initial magnetization 

We consider first the standard case where the system is 
initially unmagnetized. This implies a = \— (3,h = 2 and 
r{t) = exp(— at). From (26) the equal-time correlation is 
given by 

Cg{t, t) = e-2»*Cg(0, 0) (1 - e-2«t) {h/2a). (31) 
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For high temperatures T > Tc = I, where a > 0, 
this converges exponentiaUy to its equihbrium value 
b/2a = 1/a; the FDR approaches the limiting value 
= T(l - a)/(2 - a/a) = 1 and the system equili- 
brates as expected. Below the critical temperature, on 
the other hand, a is negative and Cg{t,t) diverges expo- 
nentially so that = 0. At criticality, finally, where 
T = 1 and o = 0, the equal-time correlator grows linearly 
as Cg{t, t) = Cg(0, 0) + 2t and the FDR has a nontrivial 
finite limit X^ = 1/2. These results can be summarized 
as follows: 




(32) 



These FDR values for the long-range Ising ferromag- 
net with zero initial magnetization are identical to those 
obtained for finite-dimensional spherical ferromagnets 
above their upper critical dimension [13, 26], as one might 
have expected on physical grounds. 



2. Non-zero initial magnetization 

For non-zero initial values of the magnetization, which 
without loss of generality we take as positive, one again 

needs to distinguish temperatures above, below and at 
the critical temperature. In the first two cases, equa- 
tion (10) tells us that m{t) decays exponentially to its 

for 



equilibrium value rricq. This value is m, 



T > Tc, while for T < Tc it is the (positive) solution of 
th(TOeq) = tanh(/3meq). Along with m, the quan- 



m, 



eq 



titles a and h also converge quickly to a^q = 1— /3(1— rrigq) 
and 6eq = 2(1 — m^q). As Oeq is just the relaxation rate 
of m{t) to meq, it is positive both above and below the 
critical temperature. From (31) the equal-time correlator 
then tends to 6eq/ (2aeq) and so the FDR approaches the 
value [32] 



^-Kq 

beq - beq/2 



1. 



(33) 



This is as expected since the system reaches equilibrium 
exponentially quickly. 

The interesting case is a quench to criticality {Tc = 1) 
with non-zero initial magnetization to(0). Here equa- 
tion (10) for the magnetization yields the asymptotic 
power-law decay 



m{t) = y/zju 



(34) 



independently of initial conditions. Also a(t) = 
tanh^(TO(<)) which at long times becomes a(t) = 3/(2t); 
as a consequence, the function r(€) from (27) scales 
asymptotically as exp[— (3/2) Int] = t~^l'^ . For the equal- 
time correlation (26), both the term including the initial 
condition Cg(0, 0) and the correction arising from the ap- 
proach of h{t) = 2[1 - m{t) tanh(TO(i))] = 2 - 3/i . . . 



to its limit value arc then subleading and one has at 
long times Cg(t, t)=2 /J dt' {t/t')-^ = t/2. The product 
a{t)Cg{t,t) thus approaches [3/{2t)] x {t/2) = 3/4 and 
the global FDR (30) tends to 



X' 



2-3/4 5 



This result can also be obtained directly from the long- 
time forms of the correlation (28) and response (29), 



Cg(t, ty,) 



-3/2 



Again, we can summarize the results: 



-3/2 



1 T < Tc 
X:^ = i 4/5 T = Tc 
1 T > Tc 



(35) 



The difi'erence between these FDR values and those for 
the unmagnetized case, eq. (32), is a clear signature of 
the difference in the underlying coarsening dynamics. For 
T < Tc it is physically obvious that the processes involved 
are very different: in the magnetized case the system 
equilibrates exponentially quickly, whereas for m = it 
ages indefinitely and equilibrium is never established. 

The result that also the FDR at criticality depends 
on whether the system is initially magnetized or not, on 
the other hand, is highly nontrivial. Indeed, one might 
have expected that the difference between the two cases 
becomes negligible at long times because the magnetiza- 
tion decays towards zero even in the initially magnetized 
scenario. Our explicit results show that this relaxation of 
m{t) does contribute significantly to the FDR, which ac- 
quires a nontrivial non-zero limiting value. The latter is 
distinct from the standard results X°° = 1/2, indicating 
that coarsening in the presence of a non-zero magneti- 
zation belongs to a different dynamical universality class 
from coarsening at m = 0. 



B. Local FDR 

Let us now check if the results obtained for global 
quantities can be reproduced from the local FDR. This 
is important because the concept of a non-equilibrium 
temperature is based on its indepencknice on the choice 
of observable, and also because numerical work has of- 
ten focussed on the simulation of local correlation and 
response functions [13, 33, 34]. 

The equal-time values of the local correlation and re- 
sponse are given in (16) and (24). Prom (17) and (23) the 
corresponding two-time quantities decay exponentially, 
so that 



Cioc(i,tw) = [l-m2(t^)]e-^ 
Rioc{t,t^) = /3[1 -tanh2(/^m(iw))]e" 



(36) 
(37) 
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where t ^ t — The FDR corresponding to the local 
correlation and response follows as 



tanh^(/3TO(tw)) 



l-m2(t^)-2m(tw)^#^ 



(38) 



and again only depends on the earlier time <w If the 
system starts in an unmagnctized state then m = at 



all times and therefore Xf, 



1 for all temperatures. 



For non-zero initial values of the magnetization and T > 
Tc, m{t) decays exponentially to zero and so Xioc(tw — *■ 
oo) — > 1. For T < Tc, m{t) also decays exponentially, 
but to a non-zero equilibrium value rngq. Nevertheless, 
because moq ~ tanh(/3moq), equation (38) implies that 
again Xioc(iw ^ oo) 1. At criticality, finally, inserting 
mQ(i) = 3/(2t) into (38) shows that also here Xioc(tw 
oo) ^ 1, though the convergence is now as a power law 
(~ 1/t^) rather than exponentially. Therefore 



for all T. 



(39) 



In summary, the limiting FDR obtained from the local 
correlation and response does not pick up any signature 
of the phase transition at Tc = 1, whether the system 
is magnetized or not. We will see in Sec. V that this is 
a somewhat pathological consequence of taking N ^ oo 
before looking at long times, and that finite-A^ correc- 
tions restore the expected correspondence between local 
and global measurements. 



The Hamiltonian of the short-range model is, by anal- 
ogy with (1), 



2d,.., 

In the interaction term the sum now runs over all n.n. 

pairs of spins; the interaction strength has been chosen 
to get the same critical temperature, = 1, in the limit 
d — > 00 as in the long-range model. The local fields are 
now given by hi = /i™* + (2(i)^^ ak instead of cq. (2), 
with the sum running over all n.n.s of i. For large d the 
field fluctuations Ahi are small, of 0(d~^/^), and so one 
can again linearize in Ahi. 

As for the long-range model, we proceed to study the 
correlation and response functions in this model. The 
general equations (7) and (8) can be used to derive the 
dynamical equations for the correlations. To arrive at 
explicit expressions, we need to analyze the correlations 
between spins and local fields. They can be expressed as 

{AhiAaj) = ^(AcTfcAcT,) = ^ ^ Ckj 



with Fourier transform 



^eiq-(r.-r.)(A/ijAa,) = (1 - u;q)Cq. 



Here 



1 



d 



COSQa 



IV. FINITE-DIMENSIONAL MODELS FOR 
LARGE d 

One would expect that the behaviour observed above 
for the long-range model should also appear in short- 
range models above their critical dimension. We there- 
fore now extend our discussion to the Ising model on a 
(i-dimcnsional hypcrcubic lattice with nearest neighbour 
(n.n.) interactions, in the limit of large d. A complication 
in this case is that there are multiple scalings of the spin 
correlation functions. For example, local correlations are 
0(1), those between n.n. spins scale as 0(l/d), and those 
between next nearest neighbours (n.n.n.) as 0{1/(P). In 
order to capture the 0(1) contribution of these corre- 
lation functions (and the corresponding responses) it is 
useful to consider the Fourier transforms 

Cq(t,iw)-E^''''^'''"'''^<^y 

I 

iiq(t,tw)=E'^''''^"''"""'^^'^(*'M- 

For example, the 2d n.n. spins with their correlations 
Ci,j{t,tw) of 0{l/d) given an overall contribution to 
Cq{t,tvj) of 0(1); the same is true for the 0{d^) n.n.n. 
spins with their 0{l/d^) correlations and so on. 



and the qa are the spatial components of the wavevector 
q. Using the large-d expansion of the free energy [35] one 
can show that equilibrium correlations involving higher 
powers of the field fluctuation, e.g. {{AhiY Aaj), have 
Fourier transforms which are suppressed by 0(1/ d) (away 
from criticality) . Following the same reasoning as for the 
long-range ferromagnet, we discard these subleading con- 
tributions. For the magnetization, this leads back to the 
expected mean-field equation of motion (10). The evolu- 
tion of the equal-time correlations follows by linearization 
and Fourier transformation of (7) as 



d_ 
dt 



Cq{t,t) = -2Cq(t,i) +2th'(m)(l - Wq)Cq(t,i) +6(i) 

(40) 

where the last term accounts for the fact that the local 
correlations Cii{t, t) = l—m'^{t) have a different equation 
of motion from the non-local ones. One can write an 
explicit expression for b{t) but this is not helpful because 
it depends itself on the Cq(i, i) which we are trying to 
find. Instead wc first solve (40) for arbitrary b{t) and then 
determine the latter such that the local correlations come 
out correctly. Generalizing the deflnitions of a, eq. (18), 
and r, eq. (27), to 

aq = l-th'(TO)(l-a;q) 



rq(t) = exp 



/ dt'aq{t'] 
Jo 
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the solution of (40) reads 



dt 



rlit') 



T,m- (41) 



The function h{t') can now be determined from the con- 
straint that J{dq) Cq{t,t) = Cii{t,t), where the short- 
hand notation / {dq) indicates the integral over q € 
[— TT, tt]** normalized by (Ztt)"^. Integrating (41) we find 
that 

l-m\t) = J {dq)rl{t)C^M+j\' j {dq) '^^b{t'). 

(42) 

For simplicity, wc focus in the following on initially un- 
correlated spins, i.e. Cq(0, 0) = 1 — to^(O). Using that 



/ 



{dq) e 



a;(l-Wq) 



" dqi 



,{x/d) cos q\ 



= •^+4d^ + - 



equation (42) then simplifies for large d to 

1 - m^t) = e-2*[l - m^{0)] + f dt' e-^^'-'\{t'). 

Jo 

For m{t) = at all times this gives the constant value 

b{t) = 2. For more general scenarios, b still converges 
to this value at long times as long as the magnetization 
decays to zero, i.e. for T >Tc. 

For the two-time correlations, linearization of (8) in 
Ahi yields the evolution equation {d/dt)Cq{t,tv,) = 
— OqC'q(i)iw) and thus 



Cq{t, ty,) 



iit) 



rq{U 



The instantanteous response remains purely local as 
in (A12), with Rq{t,t) = th'(m) to leading order in 1/d. 
The evolution of the two-time response is obtained as in 
the long-range case, by linearizing (A8) in the applied 
field. This gives simply {d / dt)Rq{t,t^) = —aqRq{t,ty,), 
which integrates to 



-^q(^? tw) 



r-q(tw) 



th'(m(iw))- 



On the basis of the correlation and response functions 
one can define a Fourier-component FDR. By analogy 
with (30) this can be simplified to 



Xq{t, ^w) — QQ^ 



TRq{t, tw) 



th'(m(<w)) 



^(i,tw) K^w) -aq(iw)Cq(iw,tw) ■ 

(43) 



The FDR for the global magnetization, which is the 
Fourier component of the spins with q = 0, is obtained 
from this expression as a special case. Let us now concen- 
trate on quenches to Tc where we found for the long-range 
model a nontrivial dependence of the results on the initial 
conditions. 



A. Zero initial magnetization 

For zero initial magnetization, m{t) = at all times 
and so aq{t) = ujq and rq{t) = exp(— a>qt). Correspond- 
ingly, the equal-time correlations (41) simplify to 

Cq(iw, tw) = exp(-2a;qiw) + — (l - e-^-"*-) . (44) 



The relevant scaling variable is clearly w = LUqt^, so we 
will focus on the limit of long times and low "frequencies" 
coq taken such that that w remains constant. Since Uq = 
q^/{2d) for small ujq, the scaling ojq ~ reflects the 

1 /2 

growing lengthscale 1/g ~ iw of the correlations as the 
system coarsens. The first term in (44), which arises from 
the decay of initial correlations, then becomes negligible 
compared to the second and we get 



Xq{t, iw) 



2 - a;q[l - exp(-2a;qtw)]/wq 



1 -I- exp(-2a;qtw) 



(45) 
(46) 



For lengthscales short compared to the time-dependent 
correlation length, where Wqfw '> 1, -'^q becomes equal 
to unity as expected because of effective equilibration on 
such short scales. For much larger lengthscales (wqtw ^ 
1), and in particular in the limit Wq ^ which gives the 
FDR for the magnetization, Xq approaches 1/2. These 
two limits can be seen in Fig. 1 and are consistent with 
the results for the long-range ferromagnet, but here we 
see that in addition one can interpolate smoothly between 
the two limits by varying the lengthscale considered. 
Similar behaviour is also found in the Ising chain [14], 
i.e. for d= 1, and in the spherical model [26]. Given that 
our calculation is based on a linearization in the local 
field fluctuations, it is not surprising that the result (46) 
can also be obtained from a Gaussian field theory [36]. 



B. Non-zero initial magnetization 

For non-zero initial magnetization and at criticality, 
the magnetization again decays according to (34). One 

then has th'{m{t^)) = 1 — 3/(2^^) to leading order 
for large times, and consequently aq(fw) = 1 — [1 — 
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FIG. 1: Dependence of the FDR on the lengthscale in the 

largc-fi short-range fcrromagnct, for initial conditions with 
zero (solid line, eq. (46)) and non-zero (dashed line, cq. (47)) 
magnetization. Shown on the a::-axis is the scaling variable 
cjqtvr ~ 9^tw, which is proportional to the squared ratio of the 
time-dependent correlation length (~ tl/^) to the lengthscale 
being probed by the chosen observable (1/g)- 



3/(2iw)](l - Wq)- This gives 



/ dt'a^it') 
Jo 



rq(i) = exp 

= g-t4-[t-(3/2)lnt](l-a;q) 
= t-(3/2)(l-Wq)g-Wqt 

and for long times, the equal-time correlations (41) follow 
as 



Cq(tw,tw) = t-3(l-.)e-2-<,*» 



lo \tw 



I N 3(1-Wq) 



-2a;q(tw-t')/,(f' 



For t„ ^ oo at fixed w = LUqi^,, the first term and the 
approach of b{t') to its limit 6 = 2 give only subleading 
corrections and one gets asymptotically 

Cq(*w,iw) = •^c(Wqfw) 

Wq 

Jo 

The Fourier-component FDR (43) therefore becomes 
X-\t,U) = 2-{l-[l-3/(2tw)](l-a;q)} 

= 2-[3/{2w) + l]J^c{w) (47) 

where we have neglected all terms that are subleading 
(~ 1/tvf) for long times. For short lengthscales, w ^ 1, 



one sees easily that Tc 1 and so Xq 1 , demonstrat- 
ing the expected equilibration. On the other hand, for 
lengthscales much larger than the time-dependent cor- 
relation length, i.e. w <C 1, one has J^c = w/2 and so 
Xq (2—3/4)"^ = 4/5. This applies in particular to the 
FDR for the magnetization (q = 0) and so is consistent 
with the result found above for the long-range model. 
Again there is a smooth lengthscale-dependence of the 
FDR that interpolates between local equilibrium and the 
nontrivial FDR, X°° = 4/5, for large lengthscales; see 
the dashed line in Fig. 1. This is rather reassuring: it 
tells us that there is nothing special about the magneti- 
zation, i.e. q = 0, even though in a magnetized system 
this is the only Fourier component that has a non-zero 
average. In physical terms, the FDR for the magnetiza- 
tion can also be observed by looking at lengthscales that 
are much smaller than the system size, as long as they are 
large compared to the time-dependent correlation length. 

We have only discussed the FDR for pure Fourier com- 
ponents here. However, if one considers more general 
correlations, say of spins across some finite range, one 
simply gets a mixture of the Fourier component FDRs. 
The resulting X(i,tw) will be a mixture of all Xq(i,tw) 
and interpolate between X(tw, iw) = 1 at equal times and 
a nontrivial asymptotic value in the limit ol 
well separated times. Following the reasoning in the one- 
dimensional case [14] , one can show that this asymptotic 
FDR is always identical to the FDR for the longest wave- 
length, i.e. Xq=Oj because the limit suppresses the 
contributions from all non-zero wavevectors. In particu- 
lar, this means that all observables that are linear in the 
spins, including local correlation and response, will give 
X°° = 1/2 for critical coarsening at zero magnetization, 
and X°° = 4/5 for the magnetized case. (Similarly to 
the long-range case discussed in the next section one can 
show that these asymptotic FDR values would only be 
observed for rather long time differences t — t^,, of order 
2'^ or larger, while for shorter time differences apparent 
equilibrium behaviour is obtained.) 



FINITE-SIZE CORRECTIONS FOR THE 
LONG-RANGE FERROMAGNET 



In the short-range model just discussed, local and 
global observables give the same limiting FDR X°°. But 
in the long-range model of Sec. Ill, this correspondence 
appears to be broken because the local correlation and 
response functions do not pick up any non-equilibrium 
effects. To analyze the origin of this discrepancy, we now 
study the 1/N corrections to Cioc (11) and iiioc (19). The 
calculations are sketched in Appendix B and we only give 
the main results here. 

Let us first consider again the unmagnetized case at 
high temperatures, T > T^. Keeping terms up to 
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0{N we find for the local correlation 



C'ioc(^) tw) — e 



1 



(l-/3)-i+/3r 



N 



-(l-/3)r 



+ 



N{1 



P) 
(48) 

wliere t = t ~ as before. The second term becomes 
dominant over the first for time differences r « T\nN, 
i.e. on a timescale that grows only logarithmically with 
the system size. (The l/7V-corrections in the square 
brackets can always be neglected; they would become 
relevant only for r ~ N, but by then the second term 
is exponentially larger than the first.) For a finite-size 
system the correction term therefore dominates the long- 
time dynamics, causing the decay rate of the correlation 
to slow from 1 to 1 — /? at t w Tin A''. The corresponding 
local response is related by FDT to the correlation; this 
is as expected because for the long times considered here 
the system is in equilibrium. A similar argument indi- 
cates that the above long-time results for T > remain 
unchanged if the initial magnetization is non-zero. 

Next we analyse the non-equilibrium dynamics at crit- 
icality {T = Tc = 1) starting from zero magnetization; 
here we would hope to retrieve from the 0{1/N) correc- 
tion terms a long-time FDR of = 1/2. One finds 

Cxoc{t,t^) = e"^ + ;^ [Cg(tw,iw) (1 - e-^) - re"^] . 

(49) 

Using the fact that Cg{t^,t^) = 2^^ -|- Cg(0,0) at criti- 
cality, the l/A/'-expansion now breaks down for ~ N, 
where the correction term becomes 0(1) rather than 
0{N~^) as the expansion assumes. For smaller values 
of the expansion remains valid, however, and the t^- 
derivative of the correlation becomes 

^aoc(i, fw) = e-^ + ^ (1 - e-^) 



N 



e ^[Cg(tw,tw) + l + r]. 



The third term can be neglected in the reliable regime 
tyi <^ N, but the second one again becomes domi- 
nant over the leading contribution for r w InN; for 
T - IniV > 1, one then has [d / dt^)C\oS,t^) = 2/N. 
The corresponding local response reads 



global one, just as the local global correspondence leads 
one to expect. If the limit ^ oo is taken before the 
long time limit, as we did in section II, then one implicitly 
discards the non-equilibrium regime. This is what leads 
to the apparent breaking of the correspondence with the 
global results. 

We now compare with the corresponding results for 
the non-equilibrium dynamics at criticality starting from 
non-zero initial magnetization. Evaluating the 
correction terms in the regime t ^ 1 where they are 
potentially relevant, we find 



2N Uw 



-3/2 



-3/2 



(50) 
(51) 



As in the unmagnetized case, the correction terms be- 
come significant for t ^InN. For larger time differences 
the local quantities become proportional to the global 
ones; as a consequence, the corrected local FDR Xi^^ 
crosses over from 1 to the global FDR X^^ = X^ = 4/5. 
Our main conclusion of this section is, therefore, that in 
any finite-size system the local-global correspondence is 
preserved. 

We note as an aside that the relevant scaling variable 
for the above crossovers in Xio^. is A/'exp(— r). Plots 
of Xioc versus iVexp(— r) would look similar to those 
in Fig. 1 for both the magnetized and unmagnetized 
cases, though one has to bear in mind that somewhat 
different quantities are being plotted: Fig. 1 refers to 
the lengthscale-dependence of the limiting FDR, whereas 
here we have a fixed short lengthscale (local correlation 
and response) and are looking at the system-size depen- 
dent crossover in time of the FDR to its limiting value. 

We have not considered the 1/iV-corrections at T < Tc 
in this section because it turns out - consistent with the 
fast equilibration when starting from non-zero magneti- 
zation - that here the l/A''-expansion breaks down al- 
ready for system ages of order InA^. 



VI. SUMMARY AND DISCUSSION 



1 



Riocit, iw) = 1 - ]^ i^sitw, iw) + 1 + T 



1 

N 



The second term becomes dominant over the first for 
T !vlnN, while the l/A'^-corrections in the first term can 
always be neglected for <C N. For t — \nN ^ 1 the 
corrected FDR is therefore Xioc{t,U) = {l/N)/{2/N) = 
1/2. To summarize, for systems that are old (tw ^ 1) but 
not yet equilibrated (tw ^ N)^ the FDR as a function 
of the time difference t crosses over from the equilibrium 
value Xioc = 1 to the non-equilibrium value X^^ = 1/2 
on a timescale t ^ \nN . In any finite-size system the 
limiting value of the local FDR therefore agrees with the 



In this paper we have solved analytically the non- 
equilibrium dynamics of the long-range Ising ferromagnet 
with Glauber dynamics, initially in the thermodynamic 
limit and then including also the leading finite-size cor- 
rections; our focus was on the correlation and response 
functions and the associated fiuctuation-dissipation ra- 
tio (FDR). We have also analysed the corresponding 
short-range model in the limit of large dimension, which 
provides useful additional insights into the lengthscale- 
dependence of the FDR. 

Our main result is that different nontrivial values of 
the limiting FDR X°° can result depending on the ini- 
tial conditions. In particular, for quenches to the critical 
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temperature we find that X°° = 1/2 in the standard sce- 
nario where the system is initially unmagnetized, while 
X°° = 4/5 if the initial magnetization is non-zero. We 
are not aware of any previous observations of a nontriv- 
ial (non-zero) value of X°° arising from initial conditions 
other than those traditionally considered; earlier stud- 
ies [28, 29] of strongly correlated initial conditions had 
always found either the standard value or X°° = 0. 

Our findings show that critical coarsening processes 
are fundamentally different depending on whether the 
system is magnetized or not, and the two cases must be 
considered as belonging to different dynamical universal- 
ity classes. One would have certainly expected such a 
distinction below T^, where in the magnetized case the 
system equilibrates rapidly. Our finding that the differ- 
ence persists even at Tc is much less obvious, seeing as 
even in the initially magnetized case the magnetization 
does decay towards zero at long times. The limiting FDR 
thus turns out to be a useful probe for distinguishing dif- 
ferent classes of non-equilibrium dynamics. Of course, 
the differences in the FDR also imply that the associated 
effective temperatures T/X°° differ between the magne- 
tized and unmagnetized cases; this emphasises that the 
same system can have different non-equilibrium effective 
temperatures depending on its initial preparation, which 
then reflect the physical differences in the ensuing non- 
equilibrium dynamics. 

One may ask whether our results are peculiar to 
the relatively simple scenarios that we have considered. 
This is not so: calculations in the spherical ferromagnet 
quenched to criticality [26] give exactly the same lim- 
iting FDR X°° = 4/5 for the magnetized case, in all 
dimensions d > A. To understand where this particular 
value comes from, one can write down a phenomcnologi- 
cal Langevin equation for the evolution of the fluctuating 
magnetization, 

^ = -m" + /i + C (52) 

where h is an external magnetic field and ^ is white noise 
with variance 0{1/N). In the case m = 0, because fluc- 
tuations around this value will be of 0{N~^^^), the non- 
linear term on the right-hand side of eq. (52) can be 
neglected. One thus recovers simple diffusive dynam- 
ics independently of n (as long as n > 1) which re- 
sults in the familiar value X°° = 1/2 for the limiting 
FDR. On the other hand, if the magnetization is initially 
non-zero then eq. (52) predicts that its average value de- 
cays as (to) ~ when there is no external field, 
whereas for /i 7^ it approaches (to) ~ /i^/". Com- 
paring with the standard scalings (to) ~ t~^^^'^^^ and 
(to) ~ h^/^ , respectively, one sees that we require n = 6 
and l/(n — 1) = B/ iyz); these two choices for n are con- 
sistent with each other only if the mean-field relation 
z = 2 — T] holds. We can then go ahead and, by lin- 
earizing eq. (52) in the small deviations of to from its 
average (to), calculate its correlation and response hmc- 
tion. This simple calculation gives for the limiting FDR 



X°° = (3n- l)/(4n-2) = {2^ + 3iyz)/{2l3 + Ai^z). In- 
serting the mean- field exponents z = 2, P — 1/2 and 
I' = 1/2, which imply n = 3, then leads to X°° = 4/5. 

Summarizing, the phenomenological Langevin equa- 
tion (52) predicts two different results for the limiting 
FDR of mean-field ferromagnets, depending on the ini- 
tial condition: 

_ r 1/2 if TO = 
g \4/5 if TO 7^0 

These are precisely the values that we found in 
our explicit calculations for the long-range and high- 
dimensional short-range models. One may then wonder 
whether this phenomenological approach can also be used 
to predict FDR values below the upper critical dimen- 
sion. The spherical model in d < 4, for example, has 
non-mean field exponents but still satisfies z = 2 — r/, so 
that the Langevin description is at least internally con- 
sistent. However, explicit calculations [26] show that the 
values of the FDR it predicts are incorrect. The Langevin 
dynamics (52) is therefore too simple to capture the full 
physics away from mean-field. It remains an open ques- 
tion whether appropriately generalized phenomenologi- 
cal descriptions can be used to rationalize FDR values in 
non- mean-field systems. 

We have also investigated in this paper the dependence 
of the limiting FDR on the observable considered. If 
the FDR and any associated effective temperature are 
physically meaningful, one would hope that e.g. local 
and global observables would lead to the same limit- 
ing FDRs. In the short-range model we showed that 
this is indeed the case, because the behaviour of both 
types of observables becomes dominated by the slowest, 
longest-wavelength Fourier modes {q « 0) in the limit 
of long times. The FDR for the Fourier modes them- 
selves showed the expected crossover between the values 
X^ — 1/2 and 4/5 for lengthscales larger than the time- 

— 1/2 

dependent correlation length {q ^ tw ) and X^ = 1 
for shorter, equilibrated, lengthscales. 

In the long-range model, we found that great care is 
needed when computing local FDRs because the limits 
N ^ 00 and t — 00 do not commute. If the 

thermodynamic limit ^ cxd is taken before the long- 
time limit (as we did in section II) one gets FDR values 
which are different from those obtained for global quan- 
tities because the non-equilibrium regime is effectively 
excluded. To find physically meaningful results, one has 
to take the long-time limit before the thermodynamic 
one; this then requires that the system size is kept large 
but finite as in section V. With this, the expected corre- 
spondence between local and global FDRs is recovered. 
These findings are not only of theoretical interest but 
also have two implications for numerical studies. First, if 
in long-range models one uses local observables to mea- 
sure non-equilibrium FDRs, simulations out to very long 
time differences will be required to obtain meaningful re- 
sults that can reveal non-equilibrium effects. Second, the 
fact that for global observables we could go directly to 
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the infinite-system size limit underlines the general mes- 
sage [14, 27] that such global quantities are much more 
robust tools for detecting FDT violations than their local 
counterparts. 

Our main finding that magnetized and unmagnetized 
coarsening processes at criticality can belong to differ- 
ent dynamical universality classes clearly deserves wider 
study in future work. Calculations for the spherical 
model [26] show that this distinction also holds true be- 
low the upper critical dimension d = 4, and provide 
explicit (and surprisingly nontrivial) predictions for the 
resulting FDR values. It would be interesting to com- 
plement this with simulation studies of Ising models in 
d = 2 and d = 3. Field-theoretic renormalization- group 
calculations [16] might also be possible for the 0{n) and 
n-vector models. For n = 1 these reduce to the Ising 
universality class, while for n ^ oo one would expect to 
recover the spherical model results; knowledge of the de- 
tailed dependence of the FDR values on the order param- 
eter dimensionality n as one interpolates between these 
two extreme values should help to round out the physical 
picture further. After the present work was completed we 
became aware that a first step in this direction has very 
recently been taken by the authors of Ref. [37] , who calcu- 
lated the FDR for the n-vector model with a magnetized 
initial state within an e-cxpansion aroimd d = 2. 

Finally, there is the possibility that there might be 
yet other initial conditions which give rise to distinct 
and nontrivial values of the limiting FDR. This does 
not seem likely, given that the obvious candidate case 
of strongly correlated but unmagnetized configurations 
gives X°° = and can thus be excluded. Nevertheless, 
a complete characterization of possible classes of non- 
equilibrium coarsening induced by different initial condi- 
tions certainly remains to be achieved. 



APPENDIX A: DYNAMICAL EQUATIONS FOR 
THE LONG-RANGE FERROMAGNET 



In this appendix we outline the explicit calculation 
of the dynamical equations for the correlations and re- 
sponses. The global correlation function is decomposed 
to leading order in the system size as 



N 



+ iN-l) 



N 



+ 0(Ar-2) 



= Cioc + C„i + 0(7V-i) 
= C^ + OiN-') 



(Al) 



using eqs. (11-13). At equal times, the local correlations 
are trivially 



Cu{t,t) =Ciocit,t) = l-m'^{t) 



(A2) 



while for the non-local correlations (7) implies to leading 

order 



-—Cni{t,t)+th'{m) 
x{{AhiAaj) + {AatAhj)) 
-^C^i{t,t) + -^th'(m)Cg(i,i). (A3) 



In the second line we have used that, for i j, 
{AhiAuj) = -^^^(AcTfeAa,) 



C\,JJJ) + 0(N-') 



N -I 



+ 



N -2 



N-1 



Cnl{t,t) 

N 



+ 0(Ar-2) 



N 



+ 0{N-'^). 



(A4) 



Defining the quantity a as in (18) and bearing in mind 
that Cg{t,t) = I — m'^{t) + Cni{t,t), we can then write 
the evolution equation (A3) for the non-local equal-time 
correlations as 



d_ 

dt 



Cni{t,t) = -2aCrA{t,t) + 2th'(m)(l - m^). 



For the two-time correlations we get similarly from (8) 
d 



dt 



Cij{t,t^) = -Cij{t,t^)+th'{m) 

x{Ahi{t)Aaj{t^))+0{N-^) 



(A5) 



which gives for the leading order of the non-local terms 
the dynamical evolution 



Ot 



-Cni{t,U)+th'{m)Cgit,U). (A6) 



For local correlations, on the other hand, the second term 
in (A5) is subleading so that to leading order 



d 

-g^Cloc{t, iw) 



"Cloc(i) ^w)- 



(A7) 



Finally we want the dynamical equations for the linear 
response functions Rij{t,t^) = 5{ai{t))/5hJ^^(t^). We 
assume that the field is applied to site j = 1 and, as 
for the correlation functions, use the appropriate scal- 
ings (19) and (20). Setting hf^^ = h and using (3) we 
can write 

^(cTi) = (th(m + Ahi + hS{hi))) - (ai). (A8) 

Here m refers to the value of the magnetization for the 
unperturbed system, while 



5{hi) 



N 



—r ^ Rki {t, 



(A9) 
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is the response function for the average vahie of the local 
field at site i. If we expand again in powers of A/ij, 
the first nontrivial term is proportional to {{Ahi)"^) and 
does not contribute to leading order. Writing (cTj) = 
m + hRii{t,tvj), equation (A8) thus becomes 

dtn 

— +h—Ra{t,t^) = th{m + h6{hi)) - m- hRa{t,ty,). 

Expanding to linear order in h then gives back at 0{h'^) 
the expected dynamical equation (10) for m, while at 
0{h) one gets 

d 

-g-^Rii{t,t^) = th'{m)S{h,} - Ra{t,U). 

In the non-local case {i ^ 1), where 6 {hi) = {Rioc + 
Rni)/N + 0(iV-2) from (A9) and (19,20), this gives to 
leading order 

— i?nl(i,iw) = th.'{m)[Rioc{t,t^)+Rn\{t,t^)]-Rnl{t,t^). 

(AlO) 

For the local case i = 1, on the other hand, the term 
proportional to S{hi) is subleading and one gets simply 

-i?ioc(t,iw). (All) 



These equations can be integrated forward in time once 
the instantaneous response is known. The latter is purely 
local, as one sees from 



Rij {t, t) 



d 



dhf- 



(th(/if * + m + A/ii) - a^) 



= 5y(th'(m + A/ii)). 
Thus, to leading order 

-Rioc(i,i) =th'(m), i?„i(i,t)=0. 



(A12) 



(A13) 



APPENDIX B: l/AT-CORRECTIONS FOR THE 
LONG-RANGE FERROMAGNET 



To calculate the l/A^-corrections to the local correla- 
tion and response we expand Cu = C\oc,o + C\oc,i/N + 
0{N-^) = Cioc + 0{N-^) and R^ = i?ioc,o + R\oc,i/N + 
0{N~'^) = R\oc + 0{N~'^); the magnetization, which en- 
ters Cioc,i, is similarly expanded as m = mo + nii/N + 
0{N^'^). The quantities Cioc,0: ^^ioc,o and mg are then 
the leading order values calculated in Sec. 11. Specif- 
ically, the global correlation and response are to lead- 
ing order Cg = Cioc,o + Cni and _Rg = i?ioc,o + ^ni as 
before; we will not try to calculate 1/A^-corrections to 
these global quantities because these would require the 
subleading corrections to the non-local terms and thus 
an accurate treatment of quantities of 0(iV~^). 



In (9) we now need to keep the quadratic term in A/ij. 
This gives for the magnetization: 



dm 
'dt 



= -m + {ti) 

= -m + th{m) + ^^^^th"{m)+0{N-^). 



From the definition of hi, ((A/ij)^) = (Ahiaj) for j ^ i, 
a quantity that we worked out in (A4). Expanding all 
quantities in the previous equation to order 0{N~^) then 
gives for the magnetization correction 



9m 1 



-ami + -Cg(i, i)th"(mo). 



This can be integrated but we will not give the explicit 
result here since it is not needed below. 

For the correlations, by expanding eq. (16) to 0{1/N) 
we arrive at 



C\oc,i{t,t) = -2ma{t)mi{t). 



(Bl) 



For t ^ tw '^e can use (A5) with i = j. Bearing in 
mind that {hi{t)Aa^{t^)) = C^i{t,t^)/N + OiN-"^), the 
0{N~^) terms give 



d_ 
dt 



C\oc,i{t,tv,) — —C\oc,i{t, <w) + th'(mo)C„i(<, fw) 



which integrates to 
C'ioc,i(tjiw) = e~^*~*''^Cioc,i(iw, i' 
+ 



(B2) 



^ dt' e-(*-*')th'(mo(t'))Cni(t', iw). 



For the response, keeping the (A/ij)^ term in the ex- 
pansion of (A12) gives 



Rij {t, t) — Si. 



th'(m) + ^^th'"(m)' 



+ 0{N-^) 



and the 0{N ^) terms show that the correction to the 
local instantaneous response is 

Rioc,i{t, t) = mith"(mo) + iCg(i, t)tW"{mo). (B3) 

For the two-time response, equation (A8) with the (A/ij)^ 
term retained becomes 

— + h—Rii{t,t^) = th{m + h6{hi)) + -{{Ahif) 

xth"{m + h5{h^)) (B4) 
- m - hRii{t, t^) + 0(Ar-2). 

To make progress we assume that the change of the vari- 
ance ((A/ij)^) caused by the field h is 0{h/N'^) rather 
than 0{h/N) and can therefore be neglected. This can 
be made plausible by looking at the instantaneous re- 
sponse: the variance ((Aai)^) is changed by an amount 



13 



of 0{h), while changes in all other covariances {A<TiA(jj) 
vanish. Thus ((A/ij)^) is indeed perturbed by a negligible 
amount 0{h/N'^), and one expects the response at later 
times to get no larger. We can therefore replace ((Aftj)^) 
by Cg{t,t)/N as before and regard it as /i-independent 
to the order in 1/N we are retaining. The 0{h) terms 
in (B4) then yield 



d_ 
dt 



th'(m) + %Mth"'(m) 



2N 

-Ra{t,t^) + 0{N 



S{hi) 



For the local case of interest here, i = I, one has S{hi) 
Rni/N + OiN-"^) from (A9). Therefore: 



d Id 

1 



th'(mo)— -Rni(i, 



— Rioc,o{t,tw) — j;^Rioc,i{t,tw) + 0{N ). 



(B5) 



The 0{N ^) terms give the desired equation of motion 
for the response correction 

d 

— iiloc,l(t,tw) = --Rloc,l(t,tw) +th'(mo)-Rni(i,iw) 
which integrates to 

Rioc,i{t,t^) = e-(*-*"')i?ioc,i(tw,iw) (B6) 
dt' e-(*-*'Hh'(mo(t'))^ni(t', iw)- 

In evaluating the above general predictions we start 
with the unmagnetized case, where mo(t) = 'rni{t) = 
at all times. Prom (Bl) the local equal- time correla- 
tion receives no correction, i.e. Cioc,i(i,i) = 0, while 
for the response Rioc,i{t,t) = —l3^Cg{t,t) from (B3). 
For T > Tc, we have seen that Cg{t,t) approaches its 
equilibrium value l/a=l/(l — /3) exponentially. With 
r{t) = exp(— ai) and (28) the global two-time correlation 
is then Cg(t,tw) = a~^exp(— ar), while its local ana- 
logue is given by (36) as Cioc,o(ijiw) = exp(— r). Thus 
the correction Cioc,i(i)iw) to the local correlation is for 
long times 



-lg-a(t'-t„) _ g-(t'-t„) 



(B7) 



(1 - /3)-ie-(i-'')" - e-" [(1 - (3)-^ + /3t] 



Combining this with the leading order term gives the 
result (48) discussed in the main text. To work 
out the corresponding correction to the response one 
notes from (B3) that the equal-time value Rioc,i{t,t) = 
—0^Cg{t,t) converges exponentially to — Also, the 
global response is Rg{t,t^) = /3exp(— ar) from (29) and 
the local one Rg{t,t^) = /?exp(— r) according to (37). 
By inserting these results into (B6) one finds 



-Rioc,i(^:*w) = /3e" 



-(1-/3)t 



/33 



1-/3 



+ /3 + /3V 



and it is easy to check that this is related by FDT to the 
correlation correction (B7) as it should be. 

Next consider the out-of-equilibrium dynamics at crit- 
icality (T = 1) starting from zero magnetization. For the 
correlation correction we use that Cg(t,t\v) — C'g(^w,^w) 
(a consequence of (28) together with r{t) = 1), while 
Cioc,o(i, iw) = exp(— r) as before. Then (B2) results in 



Cloc,l(i) ^w) 



dt' e 



C (t t ) - e"(*'"*") 



= Cg(tw,M(l-e"0 -^e-" 

and adding the leading order term given (49). The; calcu- 
lation for the response correction proceeds similarly, with 
Rg{t,ty,) = 1 from (29) and Rioc{t,t^) = exp(— r), and 
leads to 



-Rloc,l(^, t 



-e '^Cg(tw, iw) + 1 — e — re 



and hence (50). 

The final case of interest is the dynamics at criticality 
but with m(0) ^ 0. To obtain the leading contribu- 
tion at long times to the correction C\oc.i{t, tw), consider 
the integral term in (B2). The non-local correlations are 
C^i{t,ty,) = (iw/2)(iw/i)^/^-exp(-T). The second term 
makes a contribution which is at most 0(1), and expo- 
nentially suppressed for r ^ 1. The first term, on the 
other hand, contributes 



^ dt' e-(*-*')C'g(i', iw) « Cg{t, iw) 



(B8) 



where we have used that for ty, ^ 1 the factor Cg{t' , t^f) 
in the integrand varies negligibly in the region t — t' = 
0(1) that contributes significantly, and the exponential 
integrates to unity for t 3> 1. Taking into account 
that the first term in (B2) is also exponentially sup- 
pressed, equation (B8) gives the leading contribution to 
Cioc,i{t,tvj) for r » 1 and ^ 1- For the response, 
very similar arguments show that the leading contribu- 
tion to the integral term in (B6) is simply Rg{t,tvf); for 
time differences t ^ 1 the remainder of the integral and 
the first term in (B6) can be disregarded. This leads to 
the results (50) and (51) given in the main text. 

ACKNOWLEDGMENTS 

A. G. wishes to thank the University of Pompeu Fabra 
for its support. F. R. is supported by the Ministe- 
rio de Eduacacion y Ciencia in Spain (BFM2001-3525), 
STIPCO (HPRN-CT-2002-00319), and by the Distincio 
de la Generalitat de Catalunya. I. P. acknowledges finan- 
cial support from DGICYT of the Spanish Government 
and from the Distincio de la Generalitat de Catalunya. 
All authors have been supported by the SPHINX ESF 
program. 



14 



[1] S. de Groot and P. Mazur. Non-equilibnurn thermody- 
namics. North Holland (1962). 

[2] R. Kubo. Rep. Progr. Phys. 29 255 (1966); R. Kubo, M. 
Toda and N. Hashitsume. Statistical Physics 11 (2nd ed.). 

Springer- Vcrlag, Berlin (1991). 
[3] J. Casas- Vazquez and D. Jou. Rep. Prog. Phys. 66 1937 
(2003). 

[4] S. Edwards and A. Mehta. J. Phys. France 50 2489 
(1989); A. Mehta and S. Edwards. Physica A 157 1091 
(1989). 

[5] J. Kurchan. J. Phys.: Condens. Matter 12 6611 (2000); 
L. Berthier and J.L. Barrat. J. Chem. Phys. 116 6228 
(2002). 

[6] L.F. Cugliandolo, J. Kurchan, and L. Peliti. Phys. Rev. 

E, 55 3898 (1997). 
[7] L.F. Cughandolo and J. Kurchan. Phys. Rev. Lett 71 173 

(1993); L.F. Cugliandolo and J. Kurchan. J. Phys. A 27 

5749 (1994). 

[8] Th.M. Niouwonlmizcn. Phys. Rev. Lett. 80 5580 (1998). 
[9] S. Fielding and P. Sollich. Phys. Rev. Lett. 88 050603 
(2002). 

[10] A. Garriga and F. Ritort. Eur. Phys. J. B 20 105 (2001). 

[11] A. Crisanti and F. Ritort. J. Phys. A 36 R181 (2003). 

[12] L.F. Cugliandolo, in Slow Relaxation and Non- 
Eqmhbrium Dynamics in Condensed Matter, Les 
Houches Session 77. J.L. Barrat et al, eds. Springer- 
Verlag (2003). 

[13] C. Godreche and J.M. Luck. J. Phys. A 33 9141 (2000). 
[14] P. Mayer, L. Berthier, J. P. Garrahan and P. Sollich. Phys. 

Rev. E 68 016116 (2003). 
[15] P. Calabrcsc and A. Gambassi. JSTAT P07013 (2004). 
[16] P. Calabrcse and A. Gambassi. J. Phys. A 38 R133 

(2005). 

[17] R.J. Glauber. J. Math. Phys. 4 294 (1963). 



[18] A.J. Bray. J. Phys. A 22 L81 (1989). 

[19] H. Rieger, J. Kisker and M. Schreckenberg. Physica A 

210 324 (1994). 
[20] A. Prados, J.J. Brey and B. Sanchez-Rey. Europhys. Lett. 

40 13 (1997). 

[21] C. Godreche and J.M. Luck. J. Phys. A 33 1151 (2000). 
[22] E. Lippiello and M. Zannetti. Phys. Rev. E 61 3369 
(2000). 

[23] L.F. Cugliandolo, J. Kurchan and G. Parisi. Journal de 
Physique I 4 1641 (1994). 

[24] C. Godreche and J.M. Luck. J. Phys.: Condens. Matter 
14 1589 (2002). 

[25] P. Mayer and P. Sollich. J. Phys. A 37 9 (2004). 

[26] A. Annibale and P. Solhch. Preprint cond-mat/0510731. 

[27] P. Mayer, L. Berthier, J.P. Garrahan and P. Sollich. Phys. 
Rev. E 70 018102 (2004). 

[28] A. Picone and M. Henkel. J. Phys. A 35 5575 (2002). 

[29] M. Henkel and G. Schiitz. J. Phys. A 37 591 (2004). 

[30] F. Ritort and P. SoUich. Adv. Phys. 52 219 (2003). 

[31] P. SoUich, S. Fielding and P. Mayer. J. Phys.: Condens. 
Matter 14 1683 (2002). 

[32] This result applies straightforwardly for T ^ 0. When 
considering the limit T — » one has to bear in mind 
that the time required to roach the asymptotic regime 
diverges as tw ~ 2/T, so formally the limit tw ^ oo has 
to be taken before T —f for (33) to apply. 

[33] C. Chatelain. J. Phys. A 36 10739 (2003). 

[34] F. Ricci-Tersenghi. Phys. Rev. E 68 065104(R) (2003). 

[35] A. Georges and J.S. Yedidia. J. Phys. A 24 2173 (1991). 

[36] P. Calabrese and A. Gambassi. Phys. Rev. E 65 066120 
(2002). 

[37] A. A. Fedorenko and S. Trimper. Preprint cond- 
mat/0507112. 



